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Abstract 

a 

■ Calculations are presented of the in-plane far-infrared conductivity of a 

o " 

d x 2_ y 2-w&ve superconductor, incorporating elastic scattering due to impurities 



and inelastic scattering due to spin fluctuations. The impurity scattering is 



5_i 

modeled by short-range potential scattering with arbitrary phase shift, while 
scattering due to spin fluctuations is calculated within a weak-coupling Hub- 
bard model picture. The conductivity is characterized by a low-temperature 
residual Drude feature whose height and weight are controlled by impurity 
scattering, as well as a broad peak centered at 4Ao arising from clean-limit 
inelastic processes. Results are in qualitative agreement with experiment de- 
spite missing spectral weight at high energies. 
74.25.Nf, 74.72.-h, 72.10.Di, 72.10.Fk 
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I. INTRODUCTION 



Infrared studies of the classic low-temperature superconductors provided evidence for the 
existence of a superconducting energy gap as well as information on the phonon-mediated 
pairing interaction.Si In particular, at a low reduced temperature T/T c , the conductivity 
<J\(oj) of an s-wave superconductor shows an onset when uj exceeds 2 A, increasing to the 
normal-state value at frequencies several times 2 A.i Further, structure in ai(u) at 2 A + cu p 
reflects peaks at uj p in the effective phonon mediated interaction a 2 F(uj) and a detailed 
measurement of (Ti(u) can in principle!! be used to determine a 2 F{uo). Thus it was hoped 
that infrared measurements of a\{uS) for the high-temperature superconducting cuprates 
would provide similar information on both the gap and the pairing mechanism. However, 
the search for evidence of an energy gap in the cuprates has proven difficult.! In spite of a 
confluence of data on high-quality samples, there has been a wide, even divergent, view on 
the interpretation of these data. 

The solid curves in Fig. [l| show experimental results for the a-axis conductivity <J\{uj) 
of an untwinned YBa2Cu306.93 crystal .i The a-axis conductivity, which avoids the chains, 
is believed to provide a probe of the properties of the CuC>2 planes. In the normal state 
0\ [lS] exhibits a much larger spectral weight at high frequencies than would be present from 
a Drude form with a frequency independent lifetime. Various interpretations of this have 
been suggested. It has been proposed^ that this behavior reflects a two-component response 
consisting of a "free-carrier" Drude piece and a mid-infrared (MIR) contribution associated 
with "bound- carriers." Alternatively, a one-component quasiparticle description in which 
the quasiparticle scattering rate r _1 = 6max(T, uS) has also been used to fit the normal state 
data.! However, a consistent fit with constant prefactor b cannot be obtained in frequency 
ranges both above and below ~ 1000 cm -1 . 1 

Likewise, in the superconducting state, various explanations of the data have been pro- 
posed. In the two-component picture, it is arguedi'iH that the "free-carriers" condense to 
form the superfluid while the "bound- carriers" remain, giving rise to the MIR structure 
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which becomes more clearly visible below T c . In this view the mean free path of the free- 
electron component is large compared to £o so that one is in the clean limit in which it 
is difficult to see the 2 A onset of <Ji(u>). Moreover, the MIR absorption acts to obscure 
any small onset feature. An alternative explanation^ assigned the structure beginning at 
500 cm -1 to an s-wave 2 A onset, giving 2A /kT c ~ 8. The problem in this case is to 
understand the low-temperature absorption clearly visible below 2 A. 

Recently, following various experiments0 which suggest that the high-temperature 
cuprates are characterized by a non-s-wave gap, there have been several explanations for the 
infrared data based on a gap with nodes. Putikka and Hirschfeld calculated the ab plane 
conductivity for a d x 2__ y 2 state over a cylindrical Fermi surface, including impurity scattering 
with arbitrary phase shift .EH Carbotte et a/.Ei have reported oi{uj) calculations for a model 
with a gap Aocos#. They include impurity scattering within the Born approximation and 
treat the dynamics within a Migdal-Eliashberg approximation in which a Pb spectrum or 
alternatively a cut-off Lorentzian are used as models for the electron-boson spectral density. 
Basov et a/.0 have discussed the infrared conductivity in YBa2Cu30g.95 crystals which have 
been ion irradiated in terms of a possible c?-wave state and carrier localization produced 
by the ion irradiation. Sauls and co-workers have calculated both the in-plane and c-axis 
conductivities of a d x 2_ y 2 superconductor including impurity scattering with arbitrary phase 
shift but neglecting inelastic scattering.il! 

Here we propose to extend previous calculations of the microwave response for a d x 2_ y 2- 
wave superconductor to the infrared regime. Our basic strategy is the same as in the 
previous work:@ we seek a minimal model consisting of a d x 2_ y 2 BCS superconductor with 
quasiparticles whose inverse lifetimes are given by the sum of elastic impurity scattering and 
inelastic spin-fluctuation scattering 1/r = l/r in + l/r imp . The inelastic scattering rate is 
calculated from a spin-fluctuation exchange interaction which depends upon the d x 2_ y 2 gap 
and was previously found to give a reasonable fit to the temperature-dependent quasiparticle 
lifetime determined by the microwave measurements.^ For the present case, the dynamic 
frequency dependence of the lifetime will be important. While a more complete theory is 



needed, we believe that it is useful to understand to what extent this minimal theory, which 
has previously been used in the microwave regime with considerable success, can provide 
a basis for understanding the infrared data. We find that a number of qualitative features 
of the data may be understood within the current framework, but that some quantitative 
discrepancies remain which may provide clues to the physics missing in the minimal model. 

In Section II we discuss the model and lay out the calculation of g\{oj,T). Section III 
contains the results and comparisons with data. Our conclusions are given in Section IV. 



II. INFRARED CONDUCTIVITY 

Within the BCS framework, the real part of the x-axis conductivity, here labeled as 
ciiyj), is given by 

ImA«(w) 

<7i(w) = , (1) 

UJ 

with 

lmA xx (uj)=ne 2 J j^duj'TT[A(p,u + u')A(p,uj')][f(uj + Lu')-f(Lu')][v x (p)f . (2) 

Here v x (p) is the electron group velocity along the x-direction and A(p, uj) is the spectral 
weight of the Nambu propagator. That is, A(p,uj) = — Im G(p, uj)/tc, with 

nf v £r° + e p r 3 + Apr 1 

G(p, u) = — — v , (3) 

uj 2 - e% - A p 

and the r* are Pauli matrices. The tilde symbols indicate inclusion of self-energy corrections: 
uj = uj — VJo) £p = £ p + 5^3) an d Ap = A p + J2i- 111 the calculations which follow we use 
e p = —2 t(cos p x + cos Py) — n and A p = aA(T)(cosp2; — cosp^), with the parameter a chosen 
such that the maximum value of the gap on the Fermi surface is A(T). For this tight-binding 
band the group velocity mentioned above is given by v x (p) = 2tsinp x . 

The effect of impurity scattering is included by allowing the electron self-energy to include 
multiple scattering from a random site potential. In this case, the self-energy is given by@ 
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iimp 



10) 



TO 



<7o(w) 



-rc 



(4) 



c 2 -Jo(w) ' 

Here T = Tii/(7rN ), c = cottfo, and goiui) = l/(rrN ) J d 2 pTr[G(p,u})]/(2n) 2 , where is 
the impurity concentration, No is the normal phase density of states, and 5q is the scattering 
phase shift. The self-energy correction to the gap function Si vanishes for a <i-wave gap, and 
in the unitary limit, c = 0, only the E contribution remains. In this case, the quasiparticle 
relaxation rate due to the impurity scattering is 



r impl^J 



-2 ImE 



imp 
1 



(5) 



Results for 7-* versus u for several values of T are shown in Fig. 0. The two smaller values 



imp 



of T were used in Ref. |Tg to fit penetration depth measurements and are used again here to 
allow comparison with those previous results. 

In order to take into account the dynamic spin-fluctuation contributions to the quasi- 
particle lifetime, we include, in addition to the impurity scattering, the imaginary part of 
the self-energy that arises from spin-fluctuation exchange 
d 2 p' 1 



Im S sf (p,w) = J 



(2vr)2 2[/(w) - 1] 



x Jim V(p -p',u- E p/ ) [n{u - E p .) + !][!- f(E p ,)} \ r° + ^-r 3 + ^ r 1 



+ Im V(p -p',u + E p r) [n(cj + E pl ) + 1] f(E p/ ) It 



:; A p' T \ 
E n > ~ 



(6) 



Here n and / are the usual Bose and Fermi functions, and 

V(q,u) 



3 ^ 2 x (q,^) 

2 1 - Ux ((1,uj) 



(7) 



with 



d 2 p f 1 
(2tt) 2 \2 



1 + 



£ 



p+q c p 



+ A p+q A p 



E p+OL E p 



P+qJ 



lo - (E p+(i - E p ) + id 



1 

+ 4 



-p+q c p 



+ Ap+qAp 



Ep+qEp 



f(E } 



P+QJ 



f{E p 



U m + (E } 



p+q 



EL) + iS 



+ 



1 - 



g p+q g p + Ap+qAp 



/(E p+q ) + f{Ep) - 1 • 
uj - (£ p+q + E p ) + iS 



(8) 



and U a phenomenological interaction parameter. 

From Eq. (|6]) it follows that The inelastic spin-fluctuation induced lifetime of a quasipar- 
ticle of energy uj and momentum p in a superconductor at temperature T can be written 

as 

Jo 



Tin \P,U) 



r 



(27T) 



1 



i - m 



dv Im V(p — p', v)5 {uj — v — E p 
ApAp/ + e p Ep 



1 + 



w+|A p /| 



uj{uj — u) 
dv Im V(p — p', v)5 {v — uj — E p 

ApAp/ + EpEp 



n{u) + l}[l-f{uj-v)} 



1 



1 + 



uj{v — uj) 
- p', v)8 {u. 

ApAp/ + EpEp' 



n(u) + l]\f(y-u)] 



+ / dv Im V(p — p', v)5 {uj + v — E p > 
Jo 



n{v)[l-f{uj + v)] 



(9) 



uj{uj + v) 

The first and third terms represent scattering processes associated with the emission and 
absorption of spin fluctuations, while the second term arises from the recombination of two 
quasiparticles to form a pair. Above T c , A p goes to zero and Eq. (Q) reduces to the usual 
normal state expression. 

Using parameter values for the two-dimensional Hubbard model ((n) = 0.85, U = 2, 
T c = 0.2 1, A = 3T C ) which were previously used to model NMR relaxation time datajl! 
calculations of the temperature dependence of the inelastic scattering lifetime r hl 1 {p, uj = T) 
were previously reported.0 Here we examine the frequency dependence of r i ~ 1 (p, uj). Figure 
[3] shows Tj^ 1 versus uj for three different Fermi surface momenta.0 At a node, r i ~ 1 (p, uj) ~ uj 3 
at low energies crossing over to an approximate linear variation when uj > 3 A(T). This 
is similar to the temperature dependence of T|~ (p, oj = T) previously reported. At low 
energies, the usual quasiparticle-quasiparticle Coulomb scattering would vary as uj 2 . Here 
the extra power of uj reflects the low-energy variation of the d x 2_ y 2-wave density of states. 
For energies larger than several A(T), r i ^ 1 (p, uj) is expected to approach its normal state 
variation which is approximately linear for strong spin-fluctuation scattering and a nearly 



nested Fermi surface. Also apparent in Fig. |3| is that for p away from a gap node r i ~ 1 (p, to) 
varies approximately as {u — A p ) 3 at low energy, crossing over to linear variation at higher 
energy. 

The effects of both the elastic and inelastic quasiparticle lifetimes are included in the 
infrared conductivity calculations which follow by adding the electron self-energies due to 
elastic and inelastic processes 

S(p, u) = £ sf (p, u) + S imp H . (10) 

This total self-energy is then used in the evaluation of Eq. (|2|) . 

To make further progress we make certain simplifying approximations for calculational 
simplicity. ReS sf (p,cj) is not included in results presented; test calculations show that 
it varies slowly enough with frequency, temperature and momentum to be absorbed into 
an effective mass renormalization. Re S imp (u;) is retained, however, to allow for direct 
comparison with Refs. O and [15] at low frequencies. It remains to evaluate Im S sf (p, tu), 



as given in Eq. (^). As in Ref. 20, the momentum sums in Eqs. ([]) and ([5]) are performed 
by evaluating the integrands on a lattice of points covering the Brillouin zone and summing 
the results. (For a more detailed discussion see Ref. Since the momentum sum in 

Eq. (§) is nested inside that of Eq. (||), evaluation of Im S sf (p,u;) is computationally rather 
costly. Performing this evaluation for each point of the momentum lattice sums required to 
evaluate Eq. (|2|) would be prohibitively costly since a separate momentum sum is required 
for each value of u/ needed to numerically evaluate the frequency integral contained therein. 
It is therefore a practical necessity to replace ImE sf (p,w) in Eq. (|2]) by an appropriate 
interpolation which captures the most important features of the momentum frequency and 
temperature dependence of Im S sf (p, to). As discussed above, the frequency and momentum 
dependent of r i ^ n 1 (p,co') (at a temperature well below T c ) is reasonably well fit with an 
interpolation which varies as (to — A p ) 3 . One may ask whether the momentum dependence 
of this quantity has any influence on the infrared conductivity. Figure [| shows two curves 
which represent the conductivity calculated with 



S sf (p^) = -z[2T i ;(p,^)]- 1 r (11) 

for two different choices of the interpolation function [t|* (p, c^)] - • One choice is a (u— A(#)) 3 
interpolation. Another choice is to replace S sf (p, u) by its value at the gap node on the 
Fermi surface. This leads to a momentum-independent interpolation which varies as a; 3 for 
uj < 3 A and uj for uj > 3 A . As may be seen, the two interpolations produce nearly identical 
results. For simplicity of interpretation then, a momentum-independent interpolation is used 
in all of the o\ [to) evaluations which follow. 

Of course, in order to calculate (J\{lS) at other than the lowest temperatures, we must 
incorporate temperature dependence in the (r* n ) _1 interpolation. Figure [5] shows r i ^ 1 (p*, uj), 
where p* is the Fermi surface wavevector at a gap node, calculated for T equal to T c , 0.8T C , 
and 0.1T C . The low temperature curve shows the cubic-to-linear crossover at 3 Ao discussed 
above. The T = 0.08T C curve also shows a crossover at u = 3 A(T), but varies more slowly 
than u 3 at low temperatures. The solid curves in Fig. |5| show the interpolations used at the 
respective temperatures in the calculations of <J\{u). 

Before giving our results, we note that vertex corrections to the conductivity have been 
neglected. While this is justified in the impurity-dominated regime, where s-wave impu- 
rity vertex corrections to current-current correlations functions vanish identically for singlet 
states, it represents a further approximation in the region where the inelastic spin fluctua- 
tions provide the dominant scattering. We ignore these corrections in what follows in order 
to get a qualitative picture of the conductivity in the minimal model. 

III. RESULTS 

Our basic results for a\(u,T) are summarized in Figs. |6|-|8]. Here ai(u,T), normalized 
to cr = (ne 2 / m) / (2T C ) « <Ji(0,T c ), is plotted versus to for various temperatures and elastic 
impurity scattering parameters. A 2Ao/kT c ratio of 6 was chosen for these plots. Figure 
|6] shows ai(u,T) for various temperatures for a unitary impurity scattering strength T = 
0.018 T c . In the normal state at T = T c , the inelastic scattering strength is dominant, leading 
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to an enhancement of the spectral weight over the simple Drude uj~ 2 behavior. As the 
temperature decreases below T c and the d x 2_ y 2 gap opens, spectral weight is transferred into 
the uj = superfluid delta-function. A narrow "residual Drude" response due to scattering 
by nodal quasiparticles remains, becoming increasingly narrow and high as T/T c decreases 
until the quasiparticle lifetime is limited by elastic impurity scattering. Figure |7| illustrates 
the effect of increasing the elastic scattering rate T which broadens this Drude-like response. 
If the impurity scattering rate in the normal state is significantly smaller than Ao, this 
feature may be shown to persist crudely out to a frequency uj* determined by the crossover 
of the elastic and the inelastic scattering rates, and corresponding roughly to the position 
of the low temperature conductivity minimum in Figs. |6] and 0. In the clean limit we find 
uj* ~ (rT^) 1 / 4 if impurity scattering is resonant, and uj* ~ (TnT c ) 1 / 2 in the weak scattering 
limit c> 1 (where = T/[c 2 + 1]). At T = the weight in the residual Drude feature is 
largest in the unitary limit, but increases considerably more rapidly at finite temperatures 
for weak scattering, as seen in Fig. || This may be understood in terms of the slower T —>■ T 2 
crossover in the penetration depth in this caseJiH reflected in the temperature dependent 
part of the condensate depletion. We also note that the effective (experimentally observable) 
c<j — > conductivity in the weak scattering limit is of order <7i mp = ne 2 /2mT]sf, much larger 
than the "universal"i! resonant d-wave result <7oo — ne 2 /itltiAq in the clean limit. 

In the resonant impurity scattering limit c = 0, a "shoulder" arises at uj = A in the 
conductivity due to scattering of quasiparticles from the scattering resonance in the density 
of states at uj = to the peak at the gap "edge."i! Although this might in principle be 
used to identify the existence of strong elastic scattering, it is smeared by finite temperature 
effects and may thus be difficult to observe. It is seen in the T = conductivity displayed in 
Fig. |6| as a kink around A . Note there is no particular structure in the d-wave conductivity 
at uj = 2 A Q . 

In the limit p — > 0, the onset at 2 A in an s-wave superconductor is suppressed by the 
inability to conserve momentum and the onset appears at 4Ao, where inelastic scattering 
leads to a four quasiparticle final state.0 For the d x 2_ y 2 case, this effect is broadened by 
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the gapless nature of the state; however, a broad peak at 4Ao, with Ao the maximum of 
the d x 2_ y 2 gap, is clearly visible above the minimum conductivity at ~ uj* . The peak shifts 
downward and becomes larger as the quasiparticle relaxation rate becomes comparable to Ao, 
but is near 4 Ao in the clean limit, which obtains at these frequencies if Hubbard parameters 
consistent with normal state data are chosen. Figure [] shows how this broad peak shifts 
with 2 A /kT c . 

To summarize, a d x 2_ y 2 BCS model in which quasiparticle damping due to impurity 
scattering and spin- fluctuation scattering is included gives rise to (1) a low- frequency Drude- 
like feature whose width and strength depend upon the impurity concentration, and (2) a 
midinfrared maximum at u of order 4A whose strength depends upon 2A /kT c and the 
strength of the inelastic scattering. 



IV. CONCLUSIONS 



The dashed lines in Fig. |10| show a possible fit of our model calculation of the a-axis 
conductivity o\ for the normal and superconducting states of the untwinned YBCO crystal 
studied in Ref. The a-axis data was chosen for comparison to minimize the influence 
of absorption due to the YBCO chain layers. Here we have taken o§ = 20, 000 Q^ 1 cm -1 , 
unitary (c = 0) impurity scattering with T = 0.1 T c , and 2A /kT c = 6. Note that the value 
of do is consistent with an effective mass of order 2 and a quasiparticle lifetime at T c of order 
T c _1 , which is what we have used. 



Figure |10| demonstrates the qualitative similarity of the data to the theory, in particular 
the beginning of the onset of the residual Drude feature, and a conductivity peak at about 
1000 cm -1 , corresponding roughly to 4Ao- There are, however, several apparent discrep- 
ancies. Firstly, our model fails to give sufficient spectral weight at high frequencies. We 
believe that this represents a failure to adequately describe the normal state at these higher 
energies, equivalent to the discrepancies encountered when comparing the Marginal Fermi 
Liquid model to experiment .@ It may be, as recent numerical calculations suggest JHEE 
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that this extra weight arises from excitations from the lower Hubbard band to the nar- 
row quasiparticle band which forms at the upper edge of the lower Hubbard band when 
the system is doped away from half-filling. Such calculations also show a distinct peak at 
roughly J ~ 1000 cm -1 in the normal state conductivity. However, a complete explana- 
tion along these lines should provide an understanding not only for the MIR peak in the 
normal state of YBa2Cus06.6 and L^-ajBa^CuO^ but also for its absence in YBa2Cu3C>7 
and Bi2Sr2CaCu20s. It has been speculated that the latter, "spin gap" materials consist in 
their normal state of preformed pairs which form a superconducting condensate below T c .i! 
Such pair correlations might then give rise even for T > T c to the inelastic MIR peak at 
4 A found in the current mean field theory, whereas in the former class of compounds the 
coherence and pairing onset temperatures coincide, so a distinct MIR peak is seen only for 
T < T c . Finally, we note that we are not able to rule out the possibility of a separate MIR 
band of electronic excitations which do not participate in superconductivity.^ 

A second discrepancy arises when we try to compare values of impurity scattering rates 
found by fitting optical data as in Fig. ITUlwith those deduced from fits to the low-temperature 



penetration depth. In the unitary scattering limit, it was found in Ref. that values of T 
of roughly 10 _3 T C were characteristic of nominally pure high quality single crystals. As seen 
in Fig. [10|, low-frequency optical data on similar samples requires a defect scattering rate 
a hundred times larger within the same model. This might be taken to suggest that our 
treatment of the impurity t-matrix is too naive at this point. Note further that in order to 
obtain a consistent picture at both microwave and infrared frequencies, it is not sufficient to 
simply work in the Born limit, El since in general increasing c at fixed normal state scattering 
rate decreases the weight in the residual Drude peak (Fig. §). 

However, it is important to recognize that the upturn in the uiiu) data, obtained by a 
Kramers-Kronig transform of reflectance data, is extremely sensitive to the determination of 
the 100% reflectance point, a difficult task since resolution decreases at the lowest frequen- 
cies. Transmittance experiment s,il which are not subject to this difficulty, appear to show a 
downturn in ai(u) in this frequency range. This might be consistent with a true excitation 
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gap, but would also be consistent with a very narrow residual Drude feature below the op- 
tical resolution limit, since as uj — > the conductivity must rise to its DC value of roughly 
2 x lO 4 ^" 1 cm^ 1 . We therefore believe that the apparent disagreement with experiment at 
low frequencies is not to be taken seriously at this time. Measurements in the far-infrared 
crossover region will be extremely useful to settle this point. 

A clue to the phase shift due to defect scattering in these materials may be provided 
by a recent measurement by Homes et al, in which a shoulder was observed at 300 cm -1 
only at the lowest measurement temperature (6 K but not at 12 K) in a Ni-doped sample of 
optimally doped YBCO.0 This small but identifiable feature is quite close to the frequency 
A where we expect a shoulder due to unitary limit scattering from the Fermi level resonance 
to the BCS gap edge, provided A /T c ~ 3-4. 

Despite the discrepancies outlined above, the qualitative fit of this model suggests that 
(1) the superconducting state of YBCO is characterized by a gap with nodes and (2) the 
inelastic scattering at energies less than several A is suppressed in the superconducting 
state. A gap with nodes can give rise to the residual Drude weight observed for T <C T c , 
and to an increased quasiparticle density as the impurity scattering increases due to pair- 
breaking processes which are present for a non-s-wave gap. This conclusion is supported by 
the irradiation experiments of Basov et a/.0 and the Ni-doping studies of Homes et a/.,0 
where the weight in the residual Drude feature clearly increases with disorder. The way in 
which the inelastic scattering is suppressed at frequencies less than several A implies that 
the dominant dynamic interaction is electronic in nature so that its spectral weight in this 
region is reduced when the gap opens. The observed structure at ~ 1000 cm -1 is consistent 
with a predicted feature at 4 Ao arising from inelastic scattering from spin fluctuations, and 
the smooth onset of this peak is also a natural consequence of a gapless state. 

The fact that this same model provides a fit to the microwave penetration depth and a 
qualitative fit to the real part of the microwave conductivity provides further support for 
the specific picture of a d x 2_ y 2 gap with an underlying spin-fluctuation dynamic interaction. 
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APPENDIX: APPROXIMATE EVALUATION OF FREQUENCY-DEPENDENT 

CONDUCTIVITY 



Here we give an approximate form of the conductivity in the current model for a d x 2_ y 2 
state suitable for ready comparison with experimental data, requiring only a single numerical 
quadrature. The exact form of the conductivity for a flat band in the absence of vertex 
corrections is 



aj_(q = 0,uj; T) 



ne 



m 



du'W(u')S ± (u'), 



(Al) 



where W{u'\ uj) = [f(uj' — uj) — /(u/)]/a> and 



S±(lu';uj) = Im / ^ cos 2 
J 2n 



uj' 



0+ 



(A2) 



Here we have defined u± = Cj±{uj') = uj' — T,q{uj' ± i0 + ), ^o± = ^±(^') = =t sgnu/ytl^ — A|, 
as well as analogous primed quantities uj'± = Cj±(uj' — uj) and ^g ± = ^o±(^' — oj). In general, 
Eq. (|A2|) is a complicated expression to evaluate due to the dependence of S on uj and the 
consequent necessity of evaluating the renormalized frequencies self-consistently "in parallel" 
with the integral. However in a sufficiently clean system, self-consistency corrections are 
important only over an energy range of negligible interest for the infrared conductivity. We 
therefore make the replacement ui± — > uj' ± i/2r(uj'), where 1/t(uj) takes the approximate 
clean limit forms 
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in the unitary limit and 



min (o.Gv/rA^.rAo/lu; 
0.08 A (M/A ) 3 
0.7U 



T n \uj\/A 

0.08 A (M/A ) 

0.7M 



uj* < \uj\ < 3 A 
\u\ > 3 A 

\u\ < UJ* 

uj* < \uj\ < 3 A 
\u\ > 3 A 



(A3) 



(A4) 



in the Born limit. The crossover frequency uj* is determined by requiring 1/r to be contin- 



uous. 



The kernel S± may then be simplified by retaining 1/r only in the denominators of 
Eq. flA~^) , leaving 

u{M' — M) — ±(l/r - 1/t'){N' - N) uj{M - M') - |(l/r + 1/t')(N' + N) 



2^ + i(l/ r -l/r') 2 



+ 



2^ + l(l/r + l/r') 



r\2 



(A5) 



Here N(uj) is the normalized density of states. In order to compare with the full numerical 
evaluation on a square lattice, we multiply the exact flat band result for a d x 2_ y 2 state over 
a cylindrical Fermi surface with the normalized tight-binding density of states, 



N(u) 



2H K MM 

ttAo V \AoJ N t (0) 
lK ( &n.\ Nt ^ 

7r /V \\u\J N t (0) 



M < A 
\uj\ > A 



(A6) 



where K is the complete elliptic integral of the first kind. Similarly, M(ou) is approximated 



by 



M{uj) 







A 



N t (0) 



\uj\ < A 
\u\ > A 



(A7) 



The tight-binding density of states N t (uj) is given by 

1 



N t (uj) 



-K .1 



uj + fl 



(A8) 



2-kH" IV V 4t 

In Eq. (|A5| ), unprimed quantities M, iV, and 1/r are evaluated at w', whereas primed 
quantities are evaluated at uj' — uj. Substitution into Eq. (|A1|) yields an approximation which 



14 



compares well to the full numerical evaluation, with an error of only a few percent, as shown 
in Fig. [□]. 

Note that the correct high-frequency limit of the conductivity in the case of inelastic 
scattering is not given by the usual Drude expression with 1/r = 1/t{uj). For example, in 
the case of a flat band and l/r(u) = bu, the correct high-frequency limit obtained from 
Eqs. O and (|A2]) is 

(ne 2 \ b/(2u) , k , 

not o-i (ne 2 /m)[(b/uu)/(l + 6 2 )]. 
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FIG. 1. Experimental results for the real part of the a-axis conductivity in the normal, 
T = 100 K (dotted line), and superconducting, T = 20 K (solid line), states of an untwinned 
YBCO crystal! 
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FIG. 2. Quasiparticle relaxation rate r imp due to impurity scattering in the unitary limit, c = 0. 
Results are shown for T = 0.0008T C (dotted line), T = 0.018T C (dashed line), T = 0.1T C (solid line). 
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FIG. 3. Quasiparticle relaxation rate t-~ (p, u>) due to spin-fluctuation scattering. Results are 
shown as a function of frequency for T = 0.1T C . The different symbols indicate values calculated 
for p at three different points along the Fermi surface: a gap node (triangles), a gap antinode 
(circles), and a point halfway between the node and the antinode (squares). The solid line shows 
an interpolation which varies as u> for uo > 3Ao and as uj 3 for uj < 3Aq. The dashed lines vary as 
[lo — A p ] 3 with p, as above, at an antinode (long dashes) and an intermediate point (short dashes). 
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FIG. 4. Real part of the in-plane conductivity in the superconducting state for T = and 
r = 0.018T C . The two curves show results for two different choices for the self energy due to 
spin-fluctuation scattering. The solid curve incorporates the full momentum dependence reflected 
in Fig. ||. The dashed curve replaces S sf (p,u;) everywhere by its value at the gap node on the 
Fermi surface. 
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FIG. 5. Quasiparticle relaxation rate r^ 1 due to spin-fluctuation scattering. Results are shown 
as a function of frequency for T = T c (circles), T = 0.8T C (squares), and T = 0.1T C (triangles). 
The solid lines show interpolations used at each temperature for the evaluation of a\(uj). 
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FIG. 6. Real part of the in-plane conductivity for T = 0.018T C at several different temperatures. 
Results are shown for T = T c (dotted line), T = 0.8T C (short dashes), T = 0.5T C (long dashes), 
and T = (solid line). 
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FIG. 7. Real part of the in-plane conductivity in the superconducting state (T = 0). Results 
are shown for T = 0.0008T C (short dashes), T = 0.018T C (solid line), and T = 0.1T C (long dashes) 
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FIG. 8. Real part of the in-plane conductivity in the superconducting state for impurity scat- 
tering in the Born, c> 1 (dashed line), and unitary, c = (solid line), limits. Curves are shown 
for T = (dotted line is the Born limit, solid line is the unitary limit) and for T = 0.5T C (short 
dashes show the Born limit, long dashes show the unitary limit) All curves are calculated using 
r/(c 2 + 1) = 0.018T C . 
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FIG. 9. Real part of the in-plane conductivity in the superconducting state calculated for two 
different choices for the superconducting gap ratio. Results are shown for A(0) = 3T C (dotted line) 
and A(0) = 4T C (dashed line). The solid line shows the experimental result for a-axis conductivity 
in the superconducting state, T = 20 K, of an untwinned YBCO crystal.! 
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FIG. 10. Real part of the in-plane conductivity calculated for the normal, T = T c and 
T = 0.018T C (dot-dashed line), and superconducting, T = 0.1T C , states compared to experimen- 
tal results. Calculated superconducting state results are shown for T = 0.0008T C (dotted line), 
r = 0.018T C (short dashes), and T = 0.1T C (long dashes). The solid lines show experimental results 
for a-axis conductivity in the normal, T = 100 K, and superconducting, T = 20 K, states of an 
untwinned YBCO crystal.! 
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FIG. 11. A comparison of the results for the real part of the in-plane conductivity in the 
superconducting state (for T = and V = 0.018T C ) obtained from evaluation of Eq. (|T|) (dashed 
line) and Eq. ( |A1[ ) (solid line). 
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